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A  model  of  the  excimer  formation  and  decay  processes  in  electron  beam 
pumped  high  density  rare  gases  has  been  developed.  The  model  has  been 
compared  with  recent  measurements  of  the  excimer  formation  and  decay  rates 
in  Xe  and  information  about  the  radiative  lifetimes  and  dominant  collision 


mechanisms  are  obtained.  Several  important  conclusions  that  have 
bearing  on  the  vuv  laser  characteristics  of  these  systems  are  discussed. 
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SUMMARY 


The  purpose  oJ  this  program  is  to  develop  an  understanding  of  the 

molecular  kinetics  pertinent  to  new  high  efficiency,  high  power  lasers. 

We  are  studying  high  pressure  gases  that  are  initially  pumped  by 

3  -9 

intense  bursts  of  electrons  (such  as  10  amperes,  500  keV,  3  x  10  sec), 
since  these  systems  could  have  efficiencies  of  more  than  10»  and  also 
short  pulse,  high  power  capability.  This  energy,  which  is  primarily 
deposited  by  creating  atomic  ions,  rapidly  collects  in  the  lowest 
molecular  excited  state.  Laser  radiation  from  this  state  is  possible 
cither  when  the  ground  state  is  repulsive,  that  is,  when  no  bound,  ground 
molecular  state  exists  (e.g.,  the  rare  gases  and  Hg,  Cd,  and  Zn),  or 
when  the  upper  laser  state  radiates  primarily  to  higher  vibrational 
levels  of  the  ground  state. 

The  high  pressure  rare  gases  show  considerable  promise  of  providing 
highly  efficient  media  for  high  powered  laser  action  in  the  ultraviolet 
wavelength  range.  However,  since  the  excimer  radiation  is  a  continuum 
with  a  bandwidth  on  the  order  of  100  A,  the  optical  gain  is  orders  of 
magnitude  smaller  for  a  given  population  inversion  than  comparable 
discrete  transitions  and,  correspondingly,  larger  upper  state  populations 
ait  required.  Therefore,  a  large  part  of  our  research  effort  has  been 


vii 


directed  at  deteiTiining  the  conditions  for  obtaining  high  excimer 
densities  with  high  efficiency.  Various  loss  processes  have  been 
suggested  which  would  limit  laser  efficiency.  These  include  excimer- 

J 

excimer  Penning  ionization,  excimer  photoionization,  quenching  by  ground 
state  atoms,  and  absorption  by  ground  state  molecules.  A  corresponding 
effort  has  been  devoted  to  determine  the  importance  of  these  loss 
processes  and  their  effect  on  laser  operation. 

The  major  result  has  been  the  development  of  a  theoretical  model 
for  excimer  formation  and  decay  following  short  pulse  electron  beam 
excitation  of  Ar  and  Xe  at  pressures  above  1  atm.  The  model  has  been 
developed  on  the  basis  of  existing  knowledge  of  electron,  ion,  and  excited 
atom  interactions  in  rare  gases,  together  with  available  experimental 
information  on  the  behavior  of  high  pressure  rare  gases  excited  by 
charged  particles.  The  basic  framework  of  the  reaction  scheme  is 
believed  to  be  well  established  although  some  details  remain  to  be  filled 
in  and  some  of  the  rate  constants  need  to  bs  determined.  The  important 
areas  where  knowledge  is  lacking  are  pointed  out  and  the  reactions  that 
are  particularly  important  in  limiting  the  performance  as  a  laser  medium 
arc  emphasized.  The  model  has  been  used  to  calculate  excimer  formation 
and  decay  rates  in  Xe  at  different  pressures  and  the  results  are  compared 
to  recent  experimental  data  obtained  under  Febetron  excitation  conditions. 

viii 
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On  the  basis  of  our  studies  we  have  reached  the  following  con- 


lusions: 

(1)  The  excimer  is  the  upper  level  of  the  lasing  transition. 

u 

j 

(2)  The  photoionization  cross  section  of  the  E^  excimer  level 

probably  exceeds  its  stimulated  emission  cross  section.  This 
3_+ 

would  make  the  E  population  act  as  an  absorber  of  the 
u 

singlet  radiation. 

1  +  3  + 

(3)  Rapid  (electron  induced)  mixing  of  the  E  and  E  states  is 

u  u 

implied  by  the  observed  decay  behavior.  This  would  be  a 
necessary  condition  to  insure  a  population  in  the  level 
sufficient  for  lasing. 

(4)  Quenching  of  the  excimers  by  ground  state  atoms  does  not 
appear  to  be  an  important  loss  mechanism. 

(5)  Exclmer-excimer  Penning  ionization  remains  a  significant 


limitation  on  laser  efficiency 


I  INTRODUCTION 


The  excitation  of  high  density  rare  gases  by  energetic  charged 

particles,  resulting  in  radiation  in  the  vacuum  ultraviolet  region,  has 

been  of  interest  to  radiation  physicists  for  a  long  time.  In  these  high 

density  collision-dominated  systems,  the  electronic  excitation  energy  is 

quickly  and  efficiently  converted  to  excited  diatomic  molecules  known 

as  excimers.  These  excimer  states  radiate  in  a  narrow  vuv  continuum 

band  by  making  a  transition  to  a  repulsive  ground  state  that  promptly 

dissociates.  These  unique  properties  of  the  excited  high  density  rare 

gases  are  important  in  a  number  of  potential  applications.  One  of  these 

is  the  vuv  dissociation  laser,  the  feasibility  of  which  has  already  been 

1-4 

demonstrated  for  Kr  and  Xe.  A  second  promising  application  is  the 

use  of  these  systems  as  efficient  incoherent  uv  sources,  which  can  be  used 

as  laser  pumps  and  for  photochemistry  studies.  A  third  laser  oriented 

application  being  investigated  is  the  use  of  the  excimer  states  to  excite 

5 

other  molecules  by  means  of  energy  transfer  collisions.  This  appears  to 
be  an  efficient  means  of  producing  Inversions  on  transitions  in  various 
atomic  and  molecular  species.  By  this  means  it  appears  possible  to  design 
and  operate  efficient  lasers  at  many  different  wavelengths. 
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To  capitalize  on  these  applications  it  is  essential  to  have  a  clear 
physical  and  quantitative  understanding  of  the  mechanisms  of  energy 
deposition,  and  of  excimer  formation  and  relaxation  in  the  pure  high 
pressure  rare  gases.  Equally  important  is  arh  understanding  of  the  state 
structure  and  radiative  properties  of  the  rare  gas  dimer.  The  objectives 
of  this  study  have  been  to  investigate  the  structure  and  the  important 
mechanisms,  delineate  the  chain  of  reactions,  and  develop  a  model  for 
these  processes  that  serves  as  the  basis  for  a  computer  code.  These 
objectives  have  been  reached,  and  although  many  details  remain  to  be 
ascertained  and  many  reaction  rates  are  unknown,  it  is  possible  to  explain 
reasonably  all  the  observations  with  one  consistent  model.  Further,  we 
have  developed  a  computer  code  that  calculates  the  excimer  density  as  a 
function  of  time  following  a  short  excitation  pulse.  The  results  have 
been  compared  with  the  experimentally  measured  decay  rates  of  Koehler  ot  al. 
and  good  agreement  is  attained  over  an  extremely  wide  pressure  range. 


II  ANALYTICAL  MODEL 


Excimer  State  Structure 

7 

Guided  by  Mulliken's  state  structure  diagram  for  Xe,,,  we  estimated 

the  electronic  state  structure  of  Ar  from  available  scattering  data, 

<6 

spectroscopic  data,  and  theoretical  calculations.  The  resulting  curves 

8 

were  presented  in  our  semiannual  technical  report.  Since  that  time, 

+  + 

the  well  depths  of  the  Ar  and  Xe  attractive  potentials  were  measured 

6  £> 

9 

in  our  laboratory. 

In  the  case  of  Ar  *  our  measurements  confirmed  the  calculations  of 
2 

10 

Gilbert  and  Wahl  that  had  been  used  in  the  construction  of  the  Ar2 

states  reported  previously,  giving  added  confidence  in  those  results. 

For  Xc,  we  can  use  the  measured  Xe„+  well  depth,  together  with  the 

* 

calculated  ground  state  potential  curve  and  the  radiative  transition  energy 

(i.e.,  X  =  1716  ^.)  to  place  limits  on  the  Xe  excimer  potential.  We  find 

that  the  excimer  potential  must  lie  within  the  limits:  6  =  0.97  eV, 

R  =  3.39  i  and  €  =  0.75  eV,  R  =  3.12  where  e  and  R  are  the  well 
e  e  e 

den.h  and  equilibrium  internuciear  separation,  respectively.  The  well 
*  * 

depths  for  both  Xe  and  Ar  could  be  measured  directly  by  scattering 
techniques,  and  this  important  measurement  should  be  made  in  the  near 
future! 
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With  respect  to  the  laser  application,  one  of  the  most  important 
characteristics  of  the  state  structure  of  Xe2,  Kr2,  and  Ar^  is  the 
existence  of  two  close- lying  radiative  excimer  states.  These  states 
are  closely  spaced,  and  the  ground  state  is  so  steep  that  the  radiation 
appears  as  a  single  band  that  cannot  be  spectrally  resolved  into  two 
components.  It  is  expected  that  radiation  from  the  upper  singlet  state 
is  fully  allowed,  whereas  the  lower  triplet  state  will  be  longer  lived. 

The  current  status  of  experimental  and  theoretical  work  on  the  lifetimes 
is  discussed  fully  in  Section  III.  The  effect  of  spin-orbit  splitting 
on  the  excimer  potential  curves  is  discussed  in  detail  in  the  following 
subsection. 

Spin-Orbit  Sublevels  of  the  Rare  Gas  Exclmers 

% 

7 

Mulliken  provided  a  qualitative  description  of  the  Xe  excimer 

5  5 

levels  arising  from  the  5p  6s  and  5p  6p  atomic  configurations.  To  improve 

the  usefulness  of  his  analysis,  we  performed  a  semiempirical  calculation 

5 

of  the  spin-orbit  splittings  in  the  5p  6s  molecular  levels.  This  calcula- 

*  *  * 

tion  has  also  been  made  for  the  excimer  levels  of  Kr  ,  Ar  ,  and  Ne  . 

•  «  2 

♦ 

(A  similar  calculation  of  the  spin-orbit  levels  in  Hg  was  described  in 

5  5 

a  previous  report.  )  The  p  s  configuration  gives  rise  to  four  atomic 

* 

fine  structure  levels  in  each  of  the  rare  gases.  In  Ne  ,  these  are 

13  3  3 

simply  described  in  Russell-Saunders  notation  by  P. ,  P„,  P  ,  and  P  . 

1  C  l  2 

4c 

In  Xe  they  would  be  called  ($)  ,  (J)  ,  (4),,  Of)..  In  forming  the 

1  U  1  A 


4 


exeimor  molecule,  say  Xe2,  each  of  these  atomic  sublevels  produces  a 
mixture  of  T.  and  fi  components  when  projected  onto  the  internuclear  axis. 

In  the  molecule  the  atomic  angular  momentum  is  no  longer  conserved,  but 
rather  only  its  internuclear  projection.  At  large  internuclear  separation 
L  or  S  remain  coupled  and  the  projection  of  J  along  the  Internuclear  axis 
is  the  conserved  quantity  (Hund's  case  c  or  d).  As  the  internuclear 
distance  is  reduced,  the  orbital  angular  momentum  uncouples  from  the  spin 
and  couples  instead  with  the  internuclear  axis  (Hund's  case  a  or  b). 

The  states  make  a  gradual  transition  from  Hund’s  cases  (c)  and  (d)  to 
cases  (a)  and  (b) . 

3  1  „ 

Decomposing  the  P  and  P  levels  into  their  E  and  H  components 
J  1 

allows  us  to  set  up  the  matrix  Hamiltonians  to  be  diagonalized.  These 

are  shown  in  Table  I.  We  use  the  notation  of  Condon  and  Shortley,** 

where  Q  is  the  spin-orbit  splitting  parameter,  G  is  the  singlet-triplet 

splitting  parameter,  and  W  is  a  new  parameter  to  designate  the  E-IT 

splitting.  The  values  of  Q  and  G  can  be  found  by  fitting  the  atomic 

12 

energy  levels  as  reported  by  Moore.  To  calculate  the  relative  energies 
of  the  molecular  spin-orbit  levels  using  the  atomic  parameters,  we  must 
assume  that  the  spin-orbit  interaction  does  not  change  drastically  as  a 
result  of  molecular  distortion.  These  changes  are  expected  to  be  small, 


particularly  for  the  low-lying  excimer  levels,  0  ,1,0  . 

u  u  u 


5 


Figures  1  and  2  show  the  relative  energies  of  the  molecular  spin- 

*  *  5.6 

orbit  levels  of  Ar0  and  Xe^  arising  from  Ar  (3p  Is)  on  Ar  (3p  )  and 

5  6 

Xe  (5p  6s)  on  Xe  (5p  ),  respectively.  The  ordering  and  approximate 

energy  separation  of  the  states  are  shown  as  a  function  of  the  E-H 

splitting  energy.  The  E-IT  splitting  itself  is  expected  to  vary  with  R 

approximately  as  the  overlap  of  one  atom  with  the  other,  or  about 

exponentially.  As  Muliiken  suggests,  we  expect  the  E  states  to  be 

u 

lower  than  the  H  states;  hence,  the  E-0  splitting  is  taken  to  be  positive. 
Correspondingly ,  the  E  states  should  be  higher  than  the  "  states,  and 
In  that  case  the  splitting  is  negative. 

Important  conclusions  from  these  calculations  include  the  following: 

(1)  The  splitting  between  the  3 E  and  ^E  excimer  levels  should  be 

u  u 

about  the  same  for  the  various  rare  gases,  about  1600  cm  *. 

(2)  The  splitting  between  the  0  and  1  sublevels  of  the  3E+ 

u  u  u 

*  -i  *  -1 

state  should  be  quite  small  at  R^:  in  Ne2  <  1  cm  ,  Ar^  <  3  cm 

*  -1  *  -1 
Kr^  <  25  cm  ,  Xe2  <  80  era 

(3)  Potential  barriers  are  expected  as  a  result  of  spin-orbit 

+  *  ,  -1  *  -1 
mixing  in  the  0  curves  of  Ne„  (350  cm  )  and  Ar  (160  cm  ). 
u  2  2 

*  -1 
The  actual  barrier  in  Ar2  appears  to  be  smaller  (60-140  cm  , 

ref  13),  because  of  long-range  electrostatic  attraction.  These 

barriers  may  have  serious  consequences  for  the  ki  letlcs  of 

excimer  formation. 
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ORBIT  LEVEL  DIAGRAM  FOR  Ar 


FIGURE  2  SPIN  ORBIT  LEVEL  DIAGRAM  FOR  Xe. 
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(4)  No  barriers  are  predicted  for  the  0  states  of  Kr  and  Xe  . 

u  2  2 

nor  are  anv  predicted  for  the  1  and  0  excimer  levels  for 

u  u 

all  four  rare  gases. 

(5)  Finally,  we  can  estimate  the  relative  transition  strengths  of 

+  1_+ 
the  various  0  and  1  levels  by  determining  the  amount  of  Z 
u  u  u 

and  character  in  the  wave  function  for  each  state.  The 
u 

results  of  this  calculation  are  discussed  In  Section  III. 


Reaction  Model  and  Rates 


The  model  of  the  flow  of  energy  from  the  atomic  and  ionic  states, 

where  It  is  initially  deposited,  to  the  excimer  states  is  based  on  the 

excited  state  structure  of  the  rare  gas  dimers.  This  model  was  discussed 

8 

in  detail  in  our  semiannual  report  and  is  merely  outlined  here. 

The  and  3 T  exclmers  are  formed  mainly  by  3-body  collisions 

3  3  1 

of  or  atoms  with  S  ground  state  atoms.  It  can  be  showrn  that 

21  -3 

for  Xe  densities  of  10  cm  or  less,  higher  order  collision  processes 
(e.g.,  4-body  collisions)  are  not  important. 

The  energy  deposition  by  high  energy  electrons  (Including  the  effect 
of  secondaries)  in  rare  gases  produces  mainly  atomic  ions.  For  example, 
in  Ar,  78%  of  all  excited  particles  are  ions  and  22%  are  excited  atoms. 

For  this  reason,  the  electron-ion  recombination  is  a  major  source  of 
excited  atoms  when  the  excitation  density  is  high  enough.  The  dissociative 
recombination  populates  the  band  of  np  levels  (4p  in  Ar,  6p  in  Xe)  lying 


10 


.lust  above  the  s  levels  from  which  the  exclmors  are  formed.  Deacti\ation 
of  the  p  levels  to  the  s  levels  occurs  by  several  mechanisms.  At  high 


pressures,  formation  of  bound  dimers,  which  predissociate  either 
spontaneously  or  by  collision  with  ground  state  atoms,  is  a  major 
mechanism.  At  high  electron  densities,  deactivation  of  both  atomic 
and  molecular  excited  levels  by  superelastic  collisions  can  dominate  the 
relaxation  to  the  4s  levels  and  the  excimer  levels.  Our  reaction  rate 
model  is  based  on  this  simplified  picture,  which  we  believe  includes  the 
major  formation  and  relaxation  processes.  The  major  uncertainty  in  the 


model  centers  on  the  role  played  by  the  electrons.  In  the  high  electron 
density  cases  required  in  the  laser  application,  it  is  obvious  that  the 


electrons  have  a  major  effect  on  the  deactivation.  This  effect  is 
difficult  to  model  realistically  beepuse  cross  sections  of  electron 
interaction  with  excited  atoms  and  molecules  are  not  known.  We  are, 
however,  in  the  process  of  increasing  Ih ?  detail  with  which  these  inter¬ 
actions  are  treated  in  our  model.  The  electron  temperature  is  critically 
dependent  on  both  excitation  and  decxcltatlon  rates  among  the  many 


close-lying  excited  levels.  In  the  current  modeling  studies  we  have 
assumed  the  electron  temperature  to  be  constant  at  1  eV  for  reasons 
presented  in  following  sections.  .Again,  we  are  in  the  process  of  includ¬ 
ing  in  the  model  a  full  calculation  of  the  electron  temperature,  including 
energy  dependent  elastic  and  inelastic  processes. 
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The  reactions  currently  included  in  the  model  are  summarized  in 

Table  II,  together  with  our  current  best  estimates  of  the  corresponding 

reaction  rates.  These  rates  are  based  on  direct  measurements,  or  are 

estimates  based  on  related  measurements,  as  Indicated  by  the  listed 

references  and  discussed  in  our  earlier  report.  Some  rates  have  been 

modified  slightly  from  our  earlier  estimates  to  provide  better  agreement 

of  our  modeling  results  with  experiments.  A  similar  model  is  being 

14 

developed  by  George,  and  preliminary  results  of  this  model  have  been 
reported. 

Modeling  Calculations 

The  equations  for  the  ionic  and  molecular  formation,  transfer,  and 

decay  processes  summarized  in  Table  II  have  been  programmed  for  the 

computer.*  Originally,  the  nonlinear  differential  equations  were 

numerically  integrated  using  a  fourth-order  Runga-Kutta  procedure, 

16,17 

Recently,  however,  we  obtained  nn  algorithm  designed  especially 

for  the  integration  of  systems  of  "stiff"  differential  equations  (i.e., 
the  dependent  variables  are  governed  by  rates  that  differ  greatly  in 
magnitude).  This  algorithm  allows  substantially  larger  integration 
step  sizes  than  the  Runga-Kutta  method  while  retaining  the  required 

18 

accuracy.  (Treanor  has  described  a  similar  procedure  in  some  detail.  ) 
The  new  integration  package  greatly  reduces  the  computer  time  required 
for  each  case. 

*  The  computer  programming  and  calculations  have  been  very  ably  carried 
out  by  Debasich  Mukherjee. 


TABLE  II 


"V 


REACTION'S  .AND  RATES  IN  PURE  XENON  KINETIC  MODEL 


Ri’uct  i  on 

Ivo. 

Reaction 

3  **nte 

(cm  /sec ,  etc . ) 

Reference 
(see  next  page) 

l 

+ 

Xe 

+  2Xe  -*  Xe0  +  Xe 

2.5  x  10-31 

a,  b,  c,  d 

0 

X°2 

*  * 

+  o  —  Xe  +  Xo 

-7 

2  x  10 

e 

1 

No 

** 

+  2Xe  -  Xe  +  Xe 

to'31 

estimate 

0 

*  ♦ 

Xc 

* 

—  Xo  +  l»'j 

1.5  x  I 07 

estimate,  f 

7 

♦  * 

X*2 

1-  Xe  -  Xo*  +  2Xe 

to*11 

estimate 

8 

*« 

Xe2 

* 

+  e  —  Xe  *  e 

2 

io"G 

est lmate 

9 

Xe 

*  2Xe  -  Xe*  +  Xe 

2 

-32 

5  x  10 

K  i  >» 

10 

-  ?Xe  +  liv 

Variable 

1 

11 

XC2 

+  e  —  2  Xe  +  e 

-9 

10 

estimate,  ,1 

1 

12 

*  * 

Xe 

*♦  + 

+  Xe  -•  Xe  +  Xe  +  e 

-10 

5  x  10 

J.  k 

i:* 

*  * 

Xe2 

*  *  + 

+  Xe0  -  Xe^  +  2Xe  +  e 

-10 

5x10 

J.  k 

14 

Xe2 

+  Xe2  —  Xe<2  +  2Xe  *  e 

5  -x  10"10 

j.  k,  1 

17 

* 

Xe 

+•  Xe  —  Xe  +  Xe  +  e 

5  x  10”10 

J.  k 

I 
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TABLE  II  (CONCLUDED) 


REFERENCES 

E.  Iv.  McDaniel,  et  nl.,  Ion-Molecule  Reactions  ('Vi ley- Intcrsclence, 
New  York,  1070},  p.  33S. 

D.  Smith,  A.  G.  Dean,  and  I.  C.  Plumb,  J,  Phys.  B  5,  213-1  (1972), 

R,  Mahan,  J.  Chem.  Phys.  4j3,  3080  (1965). 

A.  P.  Vitols  and  H.  J.  Oskam,  Phys,  Rev.  A  (to  be  published). 

•I,  N.  Bards  ley  and  M.  A.  Biondi,  In  Advances  in  Atomic  and 
Molecular  Physics  (Academic  Press,  New  Y’ork,  1970),  Chapter  I. 

Estimate  based  on  corresponding  transitions  in  argon.  See 
",  L.  V.iese,  M.  W.  Smith,  and  B.  M.  Glennon,  Atomic  Transition 
P-obabilitie«,  Vol.  II,  NSRDS-NBS  4,  1966. 

R.  Boucique  and  P.  Mortier,  J.  Phys.  D  3,  1905  (1970), 

This  value  is  based  on  our  analysis  of  the  data  of  Koehler,  ct  al., 
Phys.  Rev.  A  (to  be  published). 

See  discussion  in  text,  Section  III. 

D,  C.  Lorents  and  R.  E.  Olson,  "Excimer  Formation  and  Decay 
Processes  in  Rare  Gases,”  Semiannual  Report  No.  1  Contract 
N00014-72-C-0457 ,  Stanford  Research  Institute,  December  1972. 

Calculation  based  on  theoretical  arguments  presented  in 
Section  IV. 

J.  B.  Gerardo  and  A.  W,  Johnson,  J.  Chem.  Phys.  59,  1738  (1973). 


Sample  calculations  of  the  time-histories  of  several  species 

concentrations  in  xenon  pumped  by  a  3-nscc,  1- Joule  electron  beam  pulse 

are  presented  in  Figures  3(a)  through  (g)  for  pressures  from  0.27  to 

6 

41.5  atm  (chosen  to  match  data  of  Koehler  et  nl.  ).  These  calculations 

* 

use  a  radiative  time  constant  for  Xe^  that  is  a  function  of  pressure 

(this  radiative  time  constant  is  discussed  fully  in  a  separate  section'. 

From  Figures  3(a)  through  (g)  we  see  that  with  the  present  model  quite 

good  agreement  can  be  achieved  botwoen  the  predicted  and  measured  con- 
* 

centration  of  Xe^,  although  it  should  be  noted  that  the  pressure 
dependent  lifetime  has  been  obtained  from  experiment  and  that  the 
experimental  data  have  been  normalized  to  the  calculations  at  the  maxi¬ 
mum  in  each  run  because  absolute  measurements  were  not  made. 

Discussion  of  the  Model 

From  a  careful  study  of  the  computer-modeled  results,  vq  find  that 

simple  analytical  relations  can  sometimes  be  used  to  predict  species 

concentrations.  For  example,  at  the  lowest  pressure  of  0,27  atm,  we 
* 

find  that  the  Xe  tlm-.  decay  is  limited  by  its  three-body  production 
« 

* 

from  Xe  (i.e.,  exclroer  formation  is  the  rate-limiting  reaction).  In 
the  following  analysis  we  also  determine  the  radiative  decay  rate  of 
Xe  .  Assuming  an  initial  reservoir  of  atoms,  [Xe  ]  ,  the  simplified 

M  V 


equations  become 


“uuuj  ,  i  ■  (rt  u*  ■*  »iUi 


drXe  ]  r  -2  P  *, 

=  -kg[Xe]  [Xe  ] 


d[Xe*] 


k9[Xe]2  [Xe*]  -  kl()[Xe*3  , 


where  kg  1b  the  three-body  association  rate  and  k  ^  is  the  radiative 
decay  rate.  These  equations  have  the  solution 

fXe  ]  =  [Xe  ]  exp  {-kgrXe]“  t]  ( 


*  ]Q  [X.]" 

[Xe  ]  = - —  (exp  {-k  [Xe]  t]  -  exp  [-k  t])  .  (4) 

kiQ  -  k9[Xe]2 

2  2 
11  kj0  »  kg[Xe]  ,  then  the  long-time  decay  rate  is  kg[Xe]  ,  and  the 

experimental  results  show  Just  this  result.  Under  these  conditions,  the 

reactions  expressed  la  equation  (2)  are  in  quasi-steady  state;  thus  the 


magnitude  of  rXe  ]  is  related  to  [Xe  ]  by 


[Xe*]  k9[Xe]" 

[Xe*]  '  k_ 


This  constant  ratio  of  populations  is  shown  in  Figure  3(a). 

Finally,  from  the  derivative  of  equation  (3),  the  maximal  [Xe  ] 
concentration  occurs  at  a  time  where 


k10  ®XP  *  (0> 
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Since  the  maximum  is  observed  to  occur  at  approximately  240  nsec  [Figure 

-i  6 

3(a)j  we  can  solve  this  implicit  equation  to  find  that  k  =  6  x  10 

Ht 

at  these  conditions.  Taking  into  account  the  fact  that  [Xe  ]  is  not 


simply  decaying  exponentially  at  the  early  times  makes  k  somewhat 
larger.  The  computer  results  shown  in  Figure  3(a)  were  generated 


assuming  k  =  9  x  10  ,  and  Indeed  the  agreement  between  the  predicted 

r 

and  measured  'Xe  i  history  is  excellent.  This  lifetime  of  110  nsec 
2 

we  believe  is  that  of  the  level. 

u 

For  high  pressures,  similar  considerations  can  be  applied  to  the 
electron  density  and  electron  temperature  at  times  after  the  atomic  ion 
concentration  has  decayed.  Then  electrons  and  molecular  ions  are 
produced  only  by  the  Penning  ionization  process  [primarily  reaction  (14), 
Table  II]  and  are  removed  by  the  recombination  reaction  (2),  Table  II. 
The  rate  equation  for  electrons  under  these  conditions  (note  that 


n  ?  rXe  ] )  is 
e  2 


e  ,  *-i2  r  -.2 

dT  ■  W  -  "2  n«  • 


Again,  a  quasi-steady  state  exists  during  the  decay,  and  the  electron 
density  is  related  to  the  excimer  concentration  by 


.(hi) 

>'J  '  a  ' 


The  constant  ratio  of  to  [Xe  ]  is  clearly  shown  in  Figures  3(e) 


through  3(g). 
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The  electron  temperature  is  governed  by  a  combination  of  several 
processes:  heeting  from  the  Penning  electrons  [reaction  (14),  Table  II], 
superelastic  collisions,  energy  loss  when  electrons  recombine  [reaction 
(2),  Table  I I ] ,  and  cooling  by  elastic  and  inelastic  collisions. 
Neglecting  the  effects  of  superelastic  and  inelastic  collisions  for  this 
discussion,  the  differential  equation  for  the  electron  temperature  is 

dTe  1  *  2 

-  =  —  [k  /XeJ  [3.3  -  T  ]  -  K  *N-n  -T  1  eV/sec  ,  (9) 

dt  n  14  2  '  e  e  e  e 

e 

where  K  is  the  rate  constant  for  elastic  cooling  and  the  electron 
e 

emitted  in  the  Penning  ionization  has  3.3  eV  energy.  If  we  again  assume 
quasi-steady  state  and  use  equation  (8)  above,  then  the  electron  tempera¬ 
ture  is 

3.3 

Te  ~  - ; - p  .  (10) 

1  +  KeN/{[Xe2](ki4k2n 

1  -5 

The  ratio  of  rates  has  the  approximate  value  K  /(k  .k„)  ~  10  (valid 

e  14  2 

for  0.2  <  T  <1  eV).  Therefore,  for  the  degree  of  excitation  ([Xe  ]/N) 

®  2 

pr^d’cted  for  Febetron  experiments,  equation  (10)  inuicates  that  the 
electron  temperature  will  vary  as  the  exclmers  decay. 

It  should  be  noted  that  the  rate  constant  for  elastic  cooling 
Increases  precipitously  for  electron  temperatures  above  1  eV,  by  as  much 
as  a  factor  of  100  at  2  eV.  Furthermore,  inelastic  collisions  begin  to 
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i 


i 

i 

f 


1 


be  significant  above  1  eV.  These  two  effects  counteract  the  rather  high 
electron  temperature  predicted  by  equation  (10),  and  should  maintain  the 
electron  temperature  at  between  1  and  1.5  eV.  At  these  electron  tempera 
tures,  electrons  play  an  important  role  in  mixing  electronic  states  of 
both  atoms  and  excimers.  This  electron  mixing  is  not  satisfactorily 
accounted  for  in  our  present  kinetic  model.  The  implications  of  the 
electron  mixing  on  exclmer  radiative  lifetime  axe  discussion  in 
Section  III. 


a 
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Til  RADIATIVE  LIFETIMES  OF  RARE-GAS  EXCIMERS 


One  of  the  fundamental  parametera  In  any  laser  system  is  the 
radiative  lifetime  of  the  upper  laser  level.  In  the  case  of  the  rare 
gas  dissociation  lasers  being  considered  here,  both  experimental  and 
analytical  efforts  thus  far  have  taken  the  simple  approach  of  considering 
one  exclmer  state  with  single  lifetime,  whereas  the  state  structure 
clearly  indicates  that  there  are  two  radiating  levels.  Measurements  of 
Intensity  decay"  *he  excimers  following  pulsed  excitation  are  difficult 
to  Interpret  on  the  basis  of  a  single  lifetime,  however,  because  results 
using  high  electron  density  pumping  differ  markedly  from  those  using  low 
electron  density  pumping,  even  at  the  same  pressure,  and  because  the 
radiative  decay  rate  in  the  high  electron  density  studies  is  pressure 
dependent.  Quenching  of  excimers  to  the  ground  Btate  has  been  advanced 
as  an  explanation  of  the  latter  effect. 

In  this  section,  we  consider  in  detail  the  contribution  of  two 
exclmer  states  to  the  observed  radiation.  These  states  have  markedly 
different  radiative  lifetimes,  and  taking  proper  account  of  the  collisions! 
mixing  of  the  states  may  furnish  an  explanation  for  the  observed  radiative 
decay  in  all  the  experiments.  After  examining  the  experimental  and 
theoretical  aspects  of  this  lifetime  problem,  we  conclude  that  the  decay 
lifetimes  observed  in  the  high  density  excitation  case  are  dominated  by 
electron  Induced  mixing  of  the  states. 


Status  of  Experiments 

Measurements  yielding  information  on  the  radiative  decay  of  rare 
gas  excimers  can  be  divided  into  two  groups.  One  group  includes  measure¬ 
ments  at  low  pressures  (<  1  atm)  with  excitation  by  very  low  current 

„  19-22 

beams  of  charged  particles  (either  electrons  or  protons).  These 

studies  indicate  the  existence  of  pressure- indpendent  decay  constants 


*  *  * 


ranging  in  magnitude  from  0.3  to  5  ^sec  for  Xe  ,  Kr  ,  Ar  ,  and  Ne  (see 

2  2  3  m 


Table  III,  p.  31).  The  second  group  includes  measurements  at  higher 


pressures  (up  to  50  atm)  obtained  with  high  current  pulsed  electron 


beam  excitation  (Febetron  706). 


1,6,23-25 


There  is,  unfortunately,  strong  disagreement  among  the  various 
experimentalists  on  how  to  interpret  and  analyze  these  high  pressure 
data.  The  problem  is  complicated  by  the  fact  that  the  excimer  decay 
results  both  from  the  radiative  decay  from  two  states  with  different 
lifetimes  and  from  an  exclmer-excimer  colllslonal  deactivatirn  (Penning 
ionization),  as  expressed  by  the  rate  equation 


d[N*]  „ 

—  '  'W  <n> 

where  rN  1  is  the  sum  of  the  concentrations  of  the  two  excimer  states, 

‘  2 

k  ,  is  the  exclner-excimer  rate,  and  k  .  is  a  weighted  combination  of  the 
14  10 

two  radiative  decay  rates  (the  treatment  of  the  decay  of  the  two  states 


as  separate  entitles  is  discussed  below).  When  the  excimer  concentration 


T 


Is  high  (at  early  times  in  the  decay  and  typically  at  high  pressure, 
since  energy  deposition  in  the  gas  increases  with  increasing  pressure), 
the  early  time  decay  of  excimer  radiation  can  be  dominated  by  the 
excimer-excimer  process.  We  believe  this  effect  is  clearly  shown  In  the 


In  subsequent  papers,  Wallace 


1.  6 

results  presented  by  Koehler  et  al.  ’ 

24 

et  al.  describe  the  decay  of  xenon  in  terms  of  early  and  late  time 

6 

constants,  and  Koehler  et  al.  find  It  necessaiy  to  fit  the  decay  of 

both  argon  and  xenon  with  the  sum  of  several  exponential  terms.  Bradley 
23 

et  al.  do  not  distinguish  two  time  constants,  but  we  believe  that 
they  also  report  an  early  time  decay  rate  that  Is  governed  by  excimer- 
excimer  processes. 

25 

On  th»  other  hand,  Gerardo  and  Johnson  have  analyzed  their 

measurements  of  the  xenon  decay  according  to  equation  (1)  and  have 

determined  the  exponential  component  of  the  decay  rato.  We  have  recon- 

« 

structed  the  xenon  data  of  Koehler  et  al.  from  their  curve  fit  parameters 
and  have  analyzed  these  data  in  terms  of  equation  (1)  to  extract  an 
exponential  decay  rate.  We  have  also  analyzed  the  few  Intensity 
histories  presented  in  references  (23)  and  (24)  and  determined  such  a 
rate.  All  these  results  are  presented  In  Figure  4.  The  various  measure¬ 
ments  are  in  reasonable  agreement  with  each  other,  and  the  exponential 
decay  rate  Is  pressure  dependent!  At  high  densities,  however,  the 
decay  rate  reaches  an  asymptotic  value.  Unfortunately,  only  the 
Livermore  results  extend  to  sufficiently  high  densities  to  show  this 

27 
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effect.  This  point  is  fundamental  to  the  interpretation  of  the  pressure 
dependence;  if  it  is  due  to  a  two-body  collisional  quenching  of  excimers 
to  the  ground  state,  the  decay  rate  will  be  linear  in  number  density  at 
all  pressures,  and  extrapolation  to  zero  pressure  gives  the  radiative 
lifetime.  The  existence  of  an  asymptote,  on  the  other  hand,  implies  a 

collisional  mixing  of  two  states  with  differing  lifetimes  (e.g., 

3  4-  1  + 

r  +  E  ).  The  pressure  dependence  of  the  decay  rate  in  this  case  ia 
u  u 

discussed  in  the  next  subsection.  Because  of  the  lmportsnce  of  a  correct 
interpretation  of  the  decay  rate,  we  hope  that  all  the  data  reported  in 
references  (23)  and  (24)  can  be  reanalyzed  to  determine  the  pure 
exponential  decay  rate  for  comparison  with  the  values  presented  in 
figure  4. 

The  various  interpretations  of  the  decay  rates,  and  their  extrapo¬ 
lation  to  zero  pressure,  have  resulted  in  estimations  of  exclmer  radiative 
* 

lifetimes  for  Xe  ranging  from  16  to  130  nsec.  The  data  of  Koehler 
6 

et  al.  indicate  that  the  decay  rates  of  argon  are  similar  to  those  of 

xenon.  These  results  should  be  compared  with  the  pressure-independent 

19  19,20 

decay  rates  of  500  nsec  for  xenon  and  2.8  and  3.7  usee  for  argon 

obtained  at  low  excitation  density. 

The  lifetime  data  of  Koeh.  or  et  al.  extends  to  pressures  as  low 
as  0,27  atm.  The  trend  of  the  data  is  toward  progressively  longer  life¬ 
times  with  decreasing  pressure  below  x  atm.  Our  kinetic  modeling  results 
show  that  at  these  low  densities  the  molecular  exclmer  production  becomes 
rate-limiting. 
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Theoretical  Predictions  of  Exoimer  Lifetimes 


The  1  sublevel  of  the  excimer  state  derives  most  of  its 
u  u 

transition  strength  to  the  ground  state  by  spin-orbit  mixing  with 

g 

the  fl  state.  This  mixing  is  represented  by 


(12 


The  state  dissociates  to  P  at  intemuclear  distance  R  *  «*>  ;  hence 
u  • 

at  very  large  R  the  coefficients  have  the  values  C  ~  1,  C  ~  C  ~  0. 

1  2  3 

At  finite  lnternuclear  distances,  Including  tho  equilibrium  distance 

R  ,  C  is  nonzero,  and  its  value  can  be  found  by  dlagonallelng  the  1 
e  J  u 

Hamiltonian  in  Table  I,  This  leads  to  the  relations 


/"icca  +  *  0  -  C3  .  0  ,  <13 

and 

2  2  2 

C1  +  C2  +  °3  =  1  *  (15 

where  l  Is  the  energy  of  the  ^  +  state. 

lu 

Assigning  a  transition  strength  of  1  to  the  pure  *P  state,  we  can 
define  the  transition  strength  of  the  3C.+  state  by 


Similarly  the  transition  strength  of  the 


V  + 


Ou 


state  would  be 


(16) 


(17) 


Calculations  of  the  relative  lifetimes  of  the  singlet  and  triplet 

excimer  states  can  now  be  made  insofar  as  we  know  the  values  of  the 

parameters  appearing  in  the  above  equations.  The  vriues  of  £  and  C  are 

known  from  atomic  parameters,  as  discussed  in  Section  II.  The  value  of 

V  depends  on  the  L-Tl  splitting  at  H  ,  which  is  unknown  but  can  be  cal- 

e 

culated  from  the  excimer  binding  energies.  We  discussed  estimates  of 

*  * 

the  binding  energies  of  Ai  and  Xe  in  Section  II;  values  for  Ne,  Ar, 

Kr,  and  Xe  are  all  believed  to  be  in  the  range  from  0,5  to  1  eV.  This 

T 

corresponds  to  E-II  splittings  of  8000  to  16000  cm  . 

The  results  of  calculations  of  the  ratio  of  triplet  ( l  )  to  singlet 

•  u 

(0*)  lifetimes  are  presented  in  Table  III  for  binding  energies  of  0,5 
and  1  eV.  The  range  of  triplet  lifetimes  is  also  given,  under  the 
assumption  that  the  (O^)  lifetime  is  about  the  same  as  the  corresponding 
atomic  state,  namely,  4  nsec.  Lifetimes  of  the  rare  gas  exc liners  that 
have  been  observed  in  various  experiments  using  low  electron  density 
pumping  are  shown  for  comparison  (we  implicitly  assume  here  that  the 
decay  rate  in  these  experiments  is  the  triplet  excimer  lifetime). 


31 


TABLE  III 


ESTIMATED  LIFETIMES  OF  TRIPLET  EXCIMERS 


Excimer 

T.  /t* 
lu  Ou 

Tlu 

If  T  +  =  4  nsec 
Ou 

Experimental 

Values 

D  =  0.5  eV 
e 

D  s  1.0 
e 

* 

- 

a 

Ne2 

1460. 

4800 

5.8  -  19  t*sec 

5.1  usee 

* 

b  c 

Ar 

2 

470 

1500 

1.9  -  6  usee 

2.8  usee,  3.7  ^sec 

* 

d  c 

Kr2 

48 

130 

0.19  -  0.52  usee 

0.3  usbc,  1.7  uaec 

* 

e 

Xe? 

24 

46 

96  -  184  nsec 

100-160  nsec, 

2 

500  neecC 

“  P.  X.  Lelchner,  Phys.  Rev.  A  8,  815  (1973). 
b 

N.  Thonnard  and  G.  S.  Hurst,  Phys.  Rev.  A  5,  1110  (1972). 
R.  Bouclque  and  P.  Mortier,  J.  Phys.  D.  3,  1905  (1970). 

d  P.  K.  Lelchner  and  R.  J.  Erlcson,  Phys.  Rev.  A  (In  press), 

e  6 

From  analysis  of  data  of  Koehler,  et  al .  at  0.27  atm. 


The  measured  values  fall  within  or  near  the  range  of  predicted  values, 

*  * 

except  for  the  1.7-usec  value  for  Kr  and  the  500-nsec  value  for  Xe  . 
The  good  agreement  in  the  other  cases  mRy  indicate  that  the  two  noted 
values  may  represent  the  radiative  lifetime  of  some  excimer  state  other 
than  the  lowest  triplet  state. 


Effect  of  Two  Lifetimes  ou  Observed  Radiative  Decay 

In  subsections  above,  we  have  described  experiments  that  indicate 

*  * 

radiative  lifetimes  for  Ar  and  Xe„  that  differ  by  as  much  as  two  orders 

2  2 

of  magnitude,  depending  on  the  experimental  conditions.  We  have  also 
indicated  the  existence  of  two  closely  spaced  excimer  states,  one  a 
triplet  and  one  a  singlet,  and  have  predicted  their  relative  radiative 
lifetimes.  In  following  paragraphs,  we  discuss  how  the  existence  of  two 
radiative  lifetimes  Influences  the  observed  radiative  decay. 

For  convenience,  we  will  orient  our  discussion  toward  Xe,  although 
it  applies  to  the  other  rare  gases  as  well.  For  Xe,  we  estimate  that  the 
fully  allowed  transition  from  the  singlet  excimer  to  the  ground  state 
has  a  radiative  lifetime  of  approximately  4  nsec,  while  the  triplet 
lifetime  should  be  In  the  range  100  to  200  nsec.  We  have  already 

19 

indicated  that  the  500  nsec  lifetime  measured  by  Bouclque  and  Mortler 
is  probably  not  the  triplet  lifetime.  On  the  other  band,  the  Febetron 
experiments  show  radiative  decay  rates  that  vary  with  pressure  in  the 


range  from  10  to  100  nsec.  In  seeking  to  explain  the  observed  results, 


we  examine  the  interplay  of  two  radiating  states  and  collisional  effects, 
Including  quenching  to  tho  ground  state  by  atoms  and  mixing  of  the  two 
states  by  atoms  and  by  electrons.  We  finally  conclude  that  electron 
mixing  dominates. 

The  rate  equations  for  decay  of  two  molecular  3tates  (neglecting 
initial  formation  processes)  are: 


dN 


dt 


Ik  +  [k 
R1  Q1 

!k  N  +  k  N  >  N 
ua  8  Ue  M 


+  k  ]  N  +  k  N  /  N 
d  a  d  ei  1 
a  e  / 


(18) 


and 

dN 
dt 

+  |k  N  +  k  N  >  N  ,  (19) 

K  “  de  *)  1 

where  N  is  population  and  the  subscripts  1,  3,  a,  and  e  refer  to  singlets, 
triplets,  ground-state  atoms,  and  electrons,  respectively.  The  sub¬ 
scripts  R,  0,  u,  and  d  on  the  rates  refer  to  radiation,  quenching  to  the 
ground  state,  and  upward  and  downward  mixing  between  the  two  states  by 
atoms  and  electrons.  The  aolutien  to  these  equations  has  the  form 


IX 


k  ]  N 


*  H  >N, 
%  ')  - 


",  -  «U.V  *  s/8*  <2°> 

and 

"a  ‘  *  ca/2*  ,  <n> 


34 


whore  the  exponential  time  constants  are 


X  =  -k  I  k  *  k  +  I  k  +  k  +  k  +  k  N  +  Ik  +  k  Ml 

1>2  ".3  L°i  °3  dn  UaJ°  LUe  deJ  *) 

±uik  -  k  +  I  k  -  k  +  k  -  k  N  +  k  •  k  K  |2 

\JR3  Ri  L°3  °i  Ua  daJ°  LUe  6J*f 

+  4|k  N  +  k  N 1  fk  N  +  k  N  1  r  (25 

LUaa  UeeJldDa  de  jJ 

We  can  now  examine  these  time  constants  under  a  variety  of  assumptions. 

At  low  N  and  N  ,  the  time  constants  simplify  to 
n  e 


*tf2-"kR1’“kR3  ’ 


which  implies  no  auenching  and  no  mixing,  and  each  state  decays  w! th 

its  own  radiative  lifetime.  If  k  »  k  ,  the  faster  component  of 

Rt  R3 

radiation  may  be  masked  by  the  excimer  formation  time,  and  the  only 

lifetime  observed  will  correspond  to  k  .  This  is  the  Justification  for 

R3 

identifying  the  100-160  nsec  decay  constant  observed  in  the  low  pressure 

Livermore  data  as  the  triplet  lifetime  for  Xe  .  Similarly,  for  Kr  , 

Z  2 

Lelchner  and  Ericson  obtained  300  nsec;  and  in  Ar,  Thonnard  and  Hurat 
found  a  lifetime  of  2,8  usee  comparable  to  Boucique  and  ttortler'a 
3.7  usee  lifetime  (Table  III). 

Now  let  us  consider  the.  effects  of  quenching  to  the  ground  state. 

If  we  neglect  k  and  kJf  the  solutions  become 
u  d 
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i.umi  n*  »i . .  wii  . . "i 


(24) 


X  =  -fk  +  k  N  ] 

1  \R1  V/' 

X  =  -fk  +  k  H  ) 

2  \R3  °3  V  • 


With  these  assumptions,  the  decay  rates  must  Increase  linearly  with 
density  at  all  densities!  In  those  experiments  where  the  pressure 


dependence  of  the  radiative  decay  is  attributed  to  quenching, 


23-25 


—1 3  3 

the  rate  k^  Is  found  to  vary  In  the  range  (0.9  -•  5)  x  10  cm  /sec. 

There  are  two  arguments  against  this  interpretation,  however.  First, 

Bouclque  and  Mortler  observed  a  500  nsec  radiative  decay  from  an  exclmer 

state  of  Xe^  that  was  Independent  of  pressure  up  to  460  torr.  This 

indicates  that  h^N  «  kR,  or  «  10  ‘  cm  /sec.  The  second,  and  more 

important,  argument  is  a  theoretical  one.  Since  there  are  no  potential 

1  + 

curve  crossings  between  the  ground  £  and  the  exclmer  states  and  since 

8 

the  separation  between  these  states  is  so  large,  collision-induced 

transitions  at  thermal  energies  should  be  completely  negligible.  These 

3  + 

states  are  analogous  to  the  N  (A  £  )  state,  which  is  quenched  by  argon 

«  u 

•17  3  —16 

at  a  rate  leas  than  10  cm  /sec,  and  by  xenon  at  less  than  2  x  10 

3 

cm  /sec  (ref.  26). 

Let  us  now  consider  the  mixing  of  singlet  and  triplet  populations 
by  atomic  collisions  (neglecting  quenching  and  mixing  by  electrons  for 
the  moment).  The  exponential  decay  rates  are 
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-J{ 


k  +  k  +  k  +  k 
"t  *3  L  d. 


].  } 

u  J  a  ) 


±j[]k  -  k  +  [k  -  k  |N  > 

\K  R,  K  V 


2  21 4 

+  4k  k  N  1* 
u  d  a  / 

a  a  / 


l ' d 


This  result  displays  the  features  of  the  experimentally  observed  decay 
rates.  In  particular  the  smaller  of  the  rates  X^  (l.e.,  the  longer, 
observed  decay  time)  Is  density  dependent,  with  the  limiting  values 


X-  -»  -k  as  N  —0  and 
2  R3 


fk  ■  k  Ik  II 

L  «3J  «.  - 


X  ^  —  k  -  f  ■ -  - — — 

2  R  k-k  +  k  +  k  N 

3  R  R  d  u  a 

L  1  3J  a  aJ 


t  .k  [_1_] .  „  [_%_] 
R„  L1  +  K  J  R>1  +  K  J 


1  c 


as  N  -«  «,  where  K  =  k  /k  Is  the  equilibrium  constant, 
a  cud 


1 

K  =  —  exp  [-(E  -  E  )/kT  ] 

c  g  13a 


\  exp  [-1000 A  ]  (30) 

ft 

Note  that  both  the  singlet  and  triplet  exclmers  decay  with  a  common  rate 
that  la  a  weighted  mixture  of  their  radiative  decay  constants. 
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We  have  replotted  the  experimentally  observed  decay  oor.stsntB  (after 

correction  for  excimer-excimer  effects)  in  Figure  5,  and  have  also  plotted 

the  solution  of  equation  (26)  for  a  set  of  parameters  that  gives  agree- 

9  -1 

rent  with  the  data.  Note  that  the  choice  of  kR  =  1.3  x  10  sec  is 

dictated  by  equations  (28)  and  (30)  together  with  the  assumed  asymptoto 

+7  -1  6 

of  th.  data  at  k  =  6  x  10  sec  .  Furthermore,  we  set  k  =  2  x  10 

3 

sec  as  discussed  above,  so  that  the  only  free  parameter  is  k  As 

V 

seen  from  Figure  5,  the  fit  of  the  curve  to  the  Febetron  experimental 

data  is  quite  appealing.  Unfortunately,  there  are  two  problems.  The 

first  Is  the  value  of  k  that  must  be  used  to  give  the  correct  asymptotic 

R1 

value.  It  corresponds  to  a  singlet  radiative  lifetime  of  0.8  nsec  and  is 

too  short  to  believe  (the  corresponding  atomic  state  has  a  lifetime  of 

4  nsec).  The  value  of  k  does  depend  on  K  ,  which  In  turn  depends  on 

l  C 

the  gas  temperature  as  indicated  by  equation  (29).  We  have  taken 

T  <  S00°K  based  on  calculations  of  the  total  energy  deposited  in  the 

Febetron  experiments.  The  temperature  would  have  to  be  substantially 

higher  for  k  to  become  a  reasonable  value.  Also,  the  value  of  k 

R1  *1 
required  to  give  the  correct  asymptotic  value  is  only  weakly  dependent 

on  k  .  For  example,  If  the  triplet  radiative  rate  is  near  the  predicted 
R3 

7-1  6-1 

value  of  10  sec  rather  than  2  x  10  sec  as  observed  by  Bouclque  and 

9 

Mortler,  then  the  value  of  k  is  reduced  only  to  1.2  x  10  .  This  Is 

Rl 


still  not  a  reasonable  value. 


The  shortcomings  of  the  explanation  outlined  ahove,  and  the  fact 
that  the  lifetime  data  generally  fall  into  categories  of  either  quite 
low  or  quite  high  electron  densities  lead  us  to  consider,  finally,  the 
mixing  of  populations  of  the  singlet  and  triplet  states  by  electrons. 

The  rates  given  by  equations  (27)  through  (30)  pertain  In  this  case  also, 
but  with  the  subscript  a  replaced  by  e,  and  we  see  that  the  lifetime  is 
a  function  of  both  electron  density  and  electron  temperature.  Thus  it 
is  possible  to  have  different  lifetimes  at  the  sane  gas  density.  This 
would  explain  not  only  the  differences  between  the  high  and  low  electron 
density  data,  but  may  also  explain  the  smaller  differences  among  the 
Febetron  results  as  noted  most  clearly  in  Figure  4.  That  is,  what 
appears  to  be  a  pressure-dependent  radiative  lifetime  could  In  fact  be 
an  electron-dependent  lifetime,  with  differences  between  various  experi¬ 
ments  reflecting  different  electron-beam  pumping  rates,  viewing  of 
different  geometric  areas  in  the  gas,  and  so  on. 

On  first  consideration,  It  would  appear  that  the  dependence  of  the 
lifetime  on  electron  density  and  temperature  would  in  fact  cause  the 
lifetime  to  vary  with  time  following  the  electron  pulse,  as  both 
electron  density  and  temperature  decay.  However,  in  the  Febetron 
experiments  the  electron  density  is  generally  high  enough  throughout  the 
decay  to  make  the  lifetime  take  lta  high-denalty  asymptotic  value  for 
the  observation  period.  Referring  to  equations  (27)  and  (28),  we  see 
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that  the  decay  rate  takes  on  its  asymptotic  value  when 

k  •  (1  +  1/K  )N  »  (k  -  k  )  .  (oi ) 

u  c  e  R  R 

e  13 

-6  3 

We  estimate  that  k  tv  10  cm  /sec  and  K  &  0.33.  If  we  take  as  an 
u  c 

e  8  -1 

upper  limit  k  =  2  x  10  sec  ,  then  we  find  that  relation  (31)  holds 

when  N  »  5  x  10  .  Our  modeling  calculations  show  that  this  inequality 

e 

holds,  for  all  but  the  lowest  pressures,  throughout  the  decay  for  the 
Febetron  experiments.  We  then  conclude  that  the  "pressure  dependence" 
observed  in  the  Febetron  experiments  is  related  to  the  electron  tempera¬ 
ture  T  . 
e 

The  electron  temperature  during  the  decay  was  discussed  in  Section 
II.  It  was  noted  that  a  high  excimer  concentrations  the  Penning 
ionization  of  exclmers  tends  to  drive  T  up  toward  3.3  eV,  but  that  the 
rapid  Increase  of  both  elastic  and  inelastic  cooling  rates  for  energies 
about  1  eV  will  counteract  the  ionization  heating  and  tend  to  keep  T 


near  1  eV.  Then  the  asymptotic  decay  rate  has  the  values 


( for  k  =  107  sec  1 ) 
Ro 


T  (eV) 
e 


k  .  ..  (sec  S 
asymptotic _ 


4.99  x  10' 
5.32  x  10 
5.44  x  107 
5.51  x  107 


5,54  x  10' 
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Thus,  the  variation  of  the  decay  rate  is  not  very  large  for  this  range 
of  electron  temperatures,  and  the  magnitude  agrees  well  with  the  values 
observed  In  the  Febetron  experiments  for  pressures  above  1  atm.  It 
remains  to  be  seen  whether  a  reasonable  variation  with  pressure  can  be 
predicted  or  if  the  decay  rate  will  be  as  nearly  constant  with  time 
throughout  the  decay  as  the  experimental  results  indicate.  We  are 
currently  modifying  our  kinetic  modeling  computer  program  so  that  we  can 
study  this  in  detail.  It  certainly  seems,  however,  that  mixing  of  the 
two  states  by  electrons  is  an  important  process  in  determining  the 


observed  exclmer  decays. 


IV  LOSS  MECHANISMS 


Several  Important  oxclmer  loss  mechanisms  exist  that  under 
certain  conditions  can  reduce  the  efficiency  for  spontaneous  radia¬ 
tion  and  for  lasing.  These  Include  colllslonal  quenching  by  ground 
state  atoms,  excimer-cxclmer  colllslonal  deactivation,  and  photo- 
absorption  by  photoionization  of  excimer  states  or  by  absorption 
from  transient  ground  state  dimers.  Our  earlier  analysis  of  the  energy 
deposition  indicated  that  an  upper  limit  to  the  efficiency  of  pro¬ 
ducing  exclmers  is  about  50$.  This  will  be  reduced  further  if  any 
of  the  above  processes  are  important.  Measurements  of  the  efficiency 
have  been  reported  ranging  from  10  to  30$.  The  photoabsorption 


processes  at?  particularly  important  for  the  laser  because  they 
reduce  the  gain  as  well  as  the  efficiency. 


Colllslonal  Losses 


Quenching  cf  the  exclmers  to  the  ground  state*  by  ground  state 

atoms  has  been  advanced  as  the  explanation  of  the  linear  pressure 

23-25 

dependence  observed  in  some  experiments,  and  in  these  cases 

-13  3 

rate  constants  on  the  order  of  10  cm  /sec  have  been  reported.  As 
discussed  in  the  previous  section,  we  believe  that  the  pressure 
dependence  is  due  instead  to  colllslonal  mixing  (mainly  by  electrons) 
of  the  singlet  and  triplet  populations.  Quenching  of  the  exclmers 


■■■'•>  H.tfk  :*c''  .V'-r -f  h 


to  the  ground  state  by  ground  state  atoms  is  extremely  slow  because 

of  the  large  separation  in  energy  between  the  ground  He  and  and  the 
3 

excited  E  at  all  lnternucloar  separations.  An  extremely  hard  colli¬ 
sion  is  required  to  distort  the  excimer  state  enough  to  permit  a  non- 
radlative  transition  to  the  ground  state.  An  example  of  a  system  that 
has  a  similar  ground  state-excited  state  separation  is  N  (a)  in  N  . 

m  m 

-17  3,  26 

In  this  case  the  quenching  rate  constant  is  less  than  10  cm  /sec. 

A  more  important  quenching  mechanism  th* .  occurs  at  high  excita¬ 
tion  densities  is  that  of  excimer -excimer  collisions.  This  is  a 
Penning  ionization  process  in  which  one  electron  is  ejected,  absorbing 
the  excess  reaction  energy,  and  one  excimer  is  deactivated.  The  rate 

constant  x’or  this  kind  of  process  has  been  measured  for  He(2^S)  + 

3  —9  3  27 

Ho(2  s)  to  be  4  x  10  cm  /esc.  Since  the  ionization  potentials  of 
#  /  3  > 

Xe2  and  of  He(2  S)  are  nearly  the  same,  we  expect  the  reaction  crosB 

sections  to  be  Blmilar,  so  that  the  rate  constant  simply  scales  with 

velocity.  On  this  basis  we  estimate  the  rate  constant  to  be  approxi- 

mately  •  8  x  10  cm  /sec;  this  should  be  a  lower  limit  since 

the  Xe  excimer  should  be  physically  larger  than  He(23s).  In  the 

-10  3  / 

calculations  reported  here,  the  value  of  k  •  5  x  10  cm  /sec  was 

14 

used,  but  the  fit  to  the  data  would  be  Improved  if  were  increased 
-9  3 

to  1.0  x  10  cm  /sec. 

Gerardo  and  Johnson  obtained  the  quantity  ki4[Xe*]Q  from  an 

44 


analysis  of  their  data. 


Using  their  50-nsec  radiative  lifetime,  they 


determined  [Xe*j  from  absolute  intensity  moasuromonts,  and  determined 
2 

-10  3  , 

k  to  be  3.5  x  10  cm  / sec.  However,  if  one  uses  the  16-nsec 

lifetime,  which  we  believe  is  correct  for  the  high  electron  density 
-9  3 

case,  then  k  ~  10  cm  /sec,  in  good  agreement  with  the  theoretical 


estimate. 


Photoionization 

28 

The  quantum  defect  method  has  been  used  for  some  years  in  the 

calculation  of  photoionization  cross  sections.  The  quantum  defect 

itBelf  enters  in  two  ways.  First,  it  describes  the  size  of  the  orbital 

occupied  by  the  electron  before  photoionization,  thus  setting  the 

general  scale  of  the  problem.  Second,  the  quantum  defect  describes 

the  phase  shift  of  the  final  state  coulomb  wave.  The  value  of  the 

dipole  matrix  element  (5  ,  ! r | *  )  can  be  quite  sensitive  to 

final  initial  n 

this  phase  shift.  This  effect  is  the  reason  for  the  small  magnitude 
of  the  threshold  photolonlzatlon  cross  sections  of  the  alkali  atoms. 

The  effect  is  washed  out  at  higher  photon  energies. 

In  estimating  the  croBS  sections  for  photoionization  of  the  rare 
gas  exclmers,  we  make  the  simple  assumption  that  the  phase  shift  of  the 
final  state  is  not  Important.  The  effects  to  be  included  are  the 
ionization  potential  (l.e.,  quantum  defect)  and  the  decrease  of  the 
dipole  matrix  element  due  to  rapid  oscillation  of  the  final  state  as 

the  impact  energy  is  Increased.  Specifically,  we  follow  the  method  of 

29  ,  . 

Bethe  and  Salpeter,  who  give  the  formula  (for  a  hydrogen  atom; 
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jto  describe  the  oBCiliator  strength,  Fn,  per  unit  of  relative  frequency, 
summed  over  all  2n2  sublevels  with  principal  quantum  number  n,  where 

2 

n2  =  ,  z  is  the  ionic  charge  after  photoioniza- 

'  '  2  J  IP  ’ 

Z  Ry 

tion,  hv  is  the  photon  energy  E,  Ry  is  the  Rydberg  unit  of  13.6  ev, 
and  IP  is  the  ionization  potential.  The  total  oscillator  strength 


per  unit  energy  is  then 


dF  $ 

n  1.96  ZRy 


where  o  is  the  relative  energy  E/lP.  The  average  oscillator  strength. 


f  jan  be  found  by  dividing  by  the  number  of  electrons  that  could  be  put 
n' 

in  the  nth  shell,  which  is  2n  (since  each  sublevel  could  be  doubly 
occupied),  to  give 


ZRy  IP  P 


Finally  we  can  calculate  the  photoionization  cross  section 


electron  as 


-17  df  ,  2. 

a  «  0.807  x  10  —  (cm  ) 


8  x  10 

-  If  " 

Z(IP/13.6)  p‘ 


(cm2) 


This  formulation  leads  to  a  calculated  threshold  photoionization 

-18  2 

cross  section  of  the  H  atom  of  8  x  10  cm  (Z  =  1,  IP=13.6,  P  =  1), 

31 

whereas  for  comparison  a  more  detailed  calculation  gives 
-18  2 

6.3  x  10  cm  .  In  particular,  we  should  note  that  for  low 
ionization  potentials  the  cross  section  Increases  Just  as  the  size  of 

4 

the  orbital,  which  is  proportional  to  IP  ,  and  falls  rapidly  as  the 

_3 

energy  is  increased  above  threshold  as  (E/lP) 

Figure  6  shows  a  comparison  of  various  experimental  measurements 
and  other  theoretical  calculations  with  the  present  theory.  In  each 
case,  the  cross  section  has  been  multiplied  by  (ip/l3,6  eV)^  and 
divided  by  the  number  of  electrons  in  the  shell  to  be  ionized,  so  that 
species  with  different  ionization  potentials  can  all  be  compared 
together.  Considering  the  simplicity  of  the  present  calculation,  the 
agreement  is  rather  good.  , 

Finally,  Figure  7  shows  our  estimates  of  the  photolonization 
cross  section  of  the  various  rare  gas  excimers,  each  evaluated  at  its 
own  emission  maximum.  The  values  are  (within  a  factor  of  2)  as  follows: 


*  „  -18 

2 

* 

-18 

*  -19 

# 

Xe^,  2  x  10 

cm  ; 

Kv 

1.3  x  10  ; 

Arg,  9  x  10  j 

2.3  x  10 

* 

and  He^,  2. 1  x 

-19 

10 

2 

cm  . 

To  arrive  at 

these  figures  we 

have 

assumed 

an  ionization  potential  of  4  eV  for  the  various  excimors.  These 
ionization  potentials  are  estimated  to  be  accurate  to  *  lo£.  A  change 
of  this  magnitude  in  the  ionization  potential  will  result  in  a  25‘t 
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FIGURE  6  COMPARISON  OF  VARIOUS  REDUCED  PHOTOIONIZATION 
CROSS  SECTIONS  WITH  PRESENT  THEORY 

He  from  reference  32,  Kr  and  Xe  reference  31,  CH4  from 
reference  33  end  34,  He  (2’S)  from  reference  38,  Ha  (E’Sfl) 
end  H2  (c3flu)  from  reference  30. 


A  •- 


change  in  the  cross  section. 


Absorption  by  Ground-State  Molecules 


Absorption  of  the  1750  A  continuum  of  xenon  by  transient  ground- 

6  . 

state  xenon  molecules  has  been  measured  by  Koehler  et  al  (measurements 

of  absorption  near  the  1470  A  atomic  resonance  line  have  also  been 
36 

reported  ).  The  results  are  anomalous  in  the  sense  that  the  tran¬ 
sient  molecular  population  is  proportional  to  the  square  of  the 
number  density  N,  and  therefore  the  absorption  coefficient  should  also 
obey  that  dependence.  In  the  pressure  range  from  0  to  30  atm,  however, 
these  measurements  show  a  linear  dependence  on  N,  namely, 

l/lQ  =  exp(  -o'  L)  (37 

-23 

where  o  ^  2.7  x  10  S.  This  linear  dependence  on  density  leads  one 
to  suspect  an  impurity.  For  example,  either  about  40  ppm  of  oxygen 
impurity  or  5  ppm  of  water  vapor  could  account  for  the  observed 
absorption  below  30  atm. 


The  molecular  xenon  absorption  is  expected  to  be  somewhat  smaller 
in  magnitude  for  pressures  less  than  30  atm  than  the  values  measured. 


We  can  show  this  approximately  by  following  the  analysis  of  Hedges 
37  38 

et  al  and  Phelps,  who  show  that  the  number  of  transient  molecules 
capable  of  absorbing  at  the  wavelength  v  (corresponding  to  an  inter- 
nuclear  distance  R)  is 


r 


\ 


;1 


x  (v)  «  N24TTR2[dv/dR)]  exp  [-V  (R)/kT] 
M  * 


(38) 


where  V^R)  is  the  repulsive  ground-state  potential,  and  dv/dR  is  the 


variation  of  absorption  frequency  with  internuclear  distance 


(-  4.5  x  1022  Hz/cm  for  xenon).  Taking  ~  0.25  eV  at  the  inter- 


nuc 


lear  distance  of  the  well  in  the  excimer  potential  (R  ~  3.3  A), 


we  find  for  T  =  300°  K  that 


.  ,  -41  2  . 

N|J(v)  ~  2  x  10  H  hv. 


(39) 


Now  the  absorption  coefficient  is 


^  X2  ^y(v) 
*(v)  -  in  \l 


5  x  10-45  N2  cm  1 


(40) 


20  -1 
Thus,  at  p  *  30  atm  (n  -  9  x  10  ),  Of  <  0.005  cm  . 


This  compares 


-1 


with  the  measured  value  of  a  b  0.025  cw 

« 

Any  of  the  absorption  mechanism*  postulated  above  (i.e,,  by 


Xe  ,  CL,  or  H  0)  can  account  for  the  observed  shift  of  wavelength  of 
2  2  2 


4,6 

the  xenon  emission  to  longer  wavelengths  with  Increasing  pressure, 
since  the  absorption  in  each  case  becomes  greater  with  decreasing 
wavelength. 

At  pressures  above  30  ata,  the  xenon  absorption  Increases  very 
rapidly.  This  is  probably  due  to  scattering  by  multi-atom  xenon 
"droplets"  that  grow  in  size  and  number  as  the  density  approaches  its 
critical  value. 


»  i 
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V  CONCLUSIONS 


A  model  of  the  excitation  of  excimers  in  high  pressure  rare 
gases  has  been  formulated  on  the  basis  of  the  electronic  structure 
proposed  for  Xe^  by  Mulllken.  This  model  includes  the  energy  deposi¬ 
tion,  electron  recombination,  exclmer  formation,  and  relaxation  pro¬ 
cesses,  and  finally,  the  radiative  properties.  Considerable  progress 
has  been  realized  in  understanding  the  observed  spontaneous  emission 
characteristics  of  the  vuv  exclmer  bands  over  a  wide  range  of  pressures 
and  excitation  densities.  The  basic  understanding  realized  in  thos j 
studies  is  directly  applicable  to  the  on-going  efforts  to  understand 
the  properties  of  these  transitions  under  lasing  conditions. 

On  the  basis  of  our  studies  we  have  reached  the  following  con¬ 
clusions. 

(1)  The  *£*  exclmer  is  the  upper  level  of  the  lasing  transition. 

(2)  The  photoionization  cross  section  of  the  exclmer  level 
probably  exceeds  its  stimulated  emission  cross  section. 

This  would  make  the  population  act  as  an  absorber  of 
the  singlet  radiation, 

(3)  Rapid  (electron  Induced)  mixing  of  the  and  states 
is  Implied  by  the  observed  decay  behavior.  This  would  be 
a  necessary  condition  to  Insure  a  population  in  the 
level  sufficient  for  lasing. 
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(4)  Excimer-excimer  Penning  ionization  remains  a  significant 
limitation  on  laser  efficiency. 

(5)  Quenching  of  the  oxcimers  by  ground  state  atoms  does  not 
appear  to  be  an  important  loss  mechanism. 

Considerable  research  remains  to  be  done  to  fully  understand  and 
to  fully  realize  the  potential  of  these  systems  as  vuv  lasers  and 
light  sources.  It  is  particularly  important  to  establish  the  elec¬ 
tronic  state  structure  of  the  excited  rare  gas  dimers,  to  determine 
rate  constants  or  cross  sections  for  many  important  collision  pro¬ 
cesses,  and  to  measure  the  radiative  transition  probabilities. 

Considerable  rate  information  is  needed  for  the  various  colli- 
slonal  relaxation  processes,  exclmer  formation  processes,  and  excimer- 
excimer  ionization  collisions.  Of  particular  importance  for  the  high 
density  excitation  cases  are  the  low  energy  electron  interactions  with 
the  exclmers  and  other  excited  states.  Electron  temperature  and 
density  measurements  of  the  plasmas  would  be  valuable. 

Radiative  lifetimes  of  the  3E  and  states  need  to  be  measured. 
The  singlet  state  lifetime  is  most  urgently  needed,  but  it  must  be 
measured  at  low  electron  density  and  low  pressure  to  avoid  collisional 
mixing  with  the  triplet  level.  Photoionization  cross  sections  of  the 
exclmer  states  should  be  measured. 

We  can  expect  that  as  more  and  better  information  is  developed 
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about  these  systems,  it  will  be  possible  to  design  efficient  and 
practical  vuv  sources,  lasers,  and  energy  transfer  media  based  on  tbs 
excimer  states.  Mainly  because  of  the  high  efficiencies  that  can  be 
attained,  it  can  be  expected  that  many  practical  applications  will  be 
realized  in  the  areas  of  new  laser  developments  and  photochemistry. 
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